Pionic BEC— BCS crossover at finite isospin chemical potential 
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Abstract 

We study the character change of the pionic condensation at finite isospin chemical potential 
\x\ by adopting the linear sigma model as a non-local interaction between quarks. At low the 
condensation is purely bosonic, then the Cooper pairing around the Fermi surface grows gradually 
as | u\ | increases. This q-q pairing is weakly coupled in comparison with the case of the q-q pairing 
that leads to color superconductivity. 
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Recent progress in computer power makes it possible to reliably simulate quantum chro- 
modynamics (QCD) at finite temperature T. As for finite density (usually parametrized 
by finite baryon chemical potential hb), however, the well known sign problem limits sim- 
ulations. Alternatively, QCD at finite isospin chemical potential fii — fi u — fid (where fi u 
and fid denoting the chemical potential of u and d quark, respectively) as well as the SU(2) 
color systems, in which the sign problem does not exist, are studied to give insights into the 
actual finite hb physics [lj]. These systems are also studied extensively in terms of effective 



models 



One of the most interesting aspects of the finite /ii systems is that they accom- 
modate pion condensation for j/^] > m n [9|, with m n denoting the mass of pions. Son and 
Stephanov [10| predicted that the pion condensed phase evolves to Cooper pairing between 
u and d (d and u) for fii > (< 0) at high but the quantitative process of the character 
change of the condensation has not been discussed. 

The BEC-BCS _ ha8 long been expected to ocenr in varions q nantn m systetns & 
[l3l | ; it was experimentally observed in ultra cold atomic gases, in which the strength of the 
interaction can be tuned artificially, only recently. At least in principle, it can occur also 
in systems governed by the strong interaction, in which the strength of the interaction 
can not be tuned artificially aside from theoretical simulations [l4[. Rather, the change in 
the environment, typically density, would lead to the crossover |15J. In symmetric nuclear 
matter, the neutron (n)-proton (p) pairing in the 3 Si + 3 Di channel that leads to bound 
deuteron formation was studied in Ref. [16] . The n-n and p-p l So pairing, that has attracted 
attention from viewpoints of both nuclear structure and neutron stars, however, does not 



reach the BEC 



17 



that the spatial extension of quar 
the mean interparticle distance 



18]. In intermediate density quark matter, the present author discussed 
c Cooper pairs in a color superconductor is comparable with 



191 ] . Later, a wide enough density region was studied [20 1 



and it was shown that the diquark pairing becomes weak at extremely high density. The 



properties of the pseudo gap phase and bosonic excitations were studied in Refs. [21 



231. 



Since the mechanism of the fermion-antifermion condensation that produces the fermion 
mass is essentially the same as the BCS pairing as recognized in Nambu and Jona-Lasinio's 



celebrated paper 



24|, the evolution of the charged pion condensation to q-q Cooper pairs 



can be analyzed in the context of the BEC-BCS crossover in terms of the spatial structure 
of the pion condensation. To this end, one must introduce a non-local interaction between 
q and q that gives momentum dependent condensations. In the present study, we adopt the 
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linear sigma model |25|], which respects chiral symmetry, as an inter-quark interaction, since 
1) the pion condensation occurs as a spontaneous symmetry breaking among three pions 
that have light but non-zero masses after the chiral symmetry breaking between the sigma 
meson and the pions, and 2) the effect of high on it has long been studied (9, [26^ 28] . 
In Ref. {29! the BEC-BCS crossover in the diquark pairing was studied in a boson-fermion 
model similar to that of the present study but the condensation is momentum independent. 
Finite jM occurs with finite /xb i n the real world; with finite T and small /xb, for example 



0.04 GeV 30J|, in heavy ion collisions and with (near) zero T and large /is, for example > 1 
GeV, in compact stars. In this sense, the present study of the system with /x B = is just the 
first step to investigate the realistic systems. However, since a signature of the BEC-BCS 
crossover in the chemical potential dependence of the condensation is measured in a lattice 



simulation for the SU(2) color system |3l( that is in a sense dual [32] to the finite /ii system, 
the spatial structure of the composite pions would be worth studying even with /j B = 0. 

When a conserved charge density J\f exists, the effective Lagrangian density is obtained 
with replacing the Hamiltonian density H by % — f/J\f, here \x denoting the correspond- 
ing chemical potential, in the partition function and performing momentum-field integra- 
tions 33]. The result for the charged pion is 

£ c g = £(7Ti -> 7Tl - /i7r 2 ,7f 2 -> 7T 2 + /J7Tl). (1) 

Since the isospin chemical potential fij corresponds to the charge chemical potential in the 
hadronic world, this form applies to the present purpose. This indicates that the role of \i\ 
corresponds to that of the angular frequency in the non-relativistic spatial rotation, that is, 
to move to a "coordinate frame" rotating in the 3 dimensional isospin space; the zero-energy 
rotational motion is a physical image of the Nambu-Goldstone mode. 
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The adopted effective Lagrangian for the quarks, sigma mesons and pions is 



£eff — £q + £m + -^couple) 

£ q = q(ifi -m q + y7°T 3 )g, 
C M = '(d^a + drf ■ d^) - U{a, t) 

+/i I (7Ti7f 2 - 7T 2 7fi) + y (vf 2 + 7r 2 ), 

U(a, 1?) = ^(a 2 + - |V + - ca, 
m 2 = A 2 / 2 - m 2 (> 0), c = /^m 2 , 

Couple = ~Gq{a + ry 5 ^ • !?)q, (2) 

where /„. and m n stand for the pion decay constant and the pion mass, respectively. Here- 
after, quantum fluctuations are indicated by primes, such as, 

q^q = {fffq) + {fffq)', 
a=(a} + a', 7T< = (n) + <■ (3) 

Since the quantum fluctuations of the quark densities and the meson fields after subtracting 
the mean field couple to each other, the normal product in C eS is understood. Note here 
that charge neutrality forced by electrons are often considered in studies of realistic /zb 



matter expected to exist in compact stars [34], |35[ . In the present study, however, charge 



neutrality is not forced since the asymmetric (/zi ^ 0) but /ib = system is an idealized one 
from the beginning. On the other hand, the charge introduced by fij is conserved among 
quarks and mesons. 

It is well known that, in the mean field level, U e g = U(a,lt) — ^-(nf + 7r|) has the 
minimum at 

<*> = ^, <-> 2 = + fl ~ <-> 2 (4) 

for l/^l > m w , assuming (n 3 ) = [9|, |26|. We take (tti) = (ir) and (7r 2 ) = without loss of 
generality. This means that the pion condensation exists in both charge sectors irrespective 
of the sign of /ij. After expanding Cm up to the quadratic terms in a' and n'^ diagonalization 
of the coupled Klein-Gordon equations for a', tt[ and 7r 2 gives the mass eigenvalues, one of 



which is zero as done in Ref. 



261 ] . But the meson mixing can not be calculated since the 



3x3 mass matrix is not regular. Thus, another approximation must be sought. Note that 



the meson mixing was calculated in another model [36] . Since the essential character of the 
massless meson propagation in the pion condensed phase is the rotational motion in the 
isospin space, we adopt a polar coordinate representation, 

7T± = -^(vrx ±m 2 ) = -^=7rexp(±z#) = -^=((tt) + vr') exp (±i6), (5) 

without expanding the angular field. This representation assures the conservation of the 
(third component of the isospin) current of the total system seen in the "rotating" frame: 

d,f = d^irfjq) + d^ x d^ 2 - t^ttO + //^(tt 2 + vr 2 ) = 0, (6) 

within the quadratic terms of the fluctuating quantum fields. In other words, the equation 
of motion of the angular field assures the current conservation. 

After confirming this point, we write down the coupled Klein-Gordon equations retaining 
the lowest order terms in each equation as 

M 0V + (2A» 2 + /i 2 )a' + 2A 2 (a)(7r>7r' = -G(qq)', 
fydV + 2A 2 (tt)V + 2\ 2 (a)(n)a' - 2// i (tt)0 = -G{qi^nq)' , 
(ir)d^e = -G(qij 5 r 2 qy, 
M 0V 3 + ^ = -G(gz 7 5 r 3 g)'. (7) 

Here we make one additional approximation to handle the set of equations: We ignore 
—2ni(7r)6 in the second equation that corresponds to the Coriolis coupling. Its influence 
will be checked later. The obtained set contains 1) the a-n mixing (the first and second 
equations), and 2) the rotational massless field (the third equation) due to the existence of 
the pion condensation (ir). 

The equation of motion of the quark propagator 

G%(x-x') = -i(tyT q i(x)^(x')\6), (8) 

where i, j and a, (3 represent isospin and Dirac indices, respectively, and |0) is the pion 
condensed ground state, is given by 

(i$ -m q + ^-1°t 3 )G(x - x') 

2 (9) 
= 5\x - x') - iG(0\T(a{x) + ij 5 l> ■ ~^(x))q(x)q(x')\0). 
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After sorting the mean field terms in 

a + 27 s • ~ (a) + ry 5 Ti(7r) + a' + i^ij^' + t 2 (ty)9 + ^71-3) (10) 

to the left-hand side, we substitute Eq.(|7D inverted by diagonalizing the meson mixing to 
Eq.flHD- Then we perform a one-body reduction (the Wick decomposition) such as 

(0\Tq(y)q(y)q(x)q(x')\0) -> (0\Tq(x)q(y)\0)(0\Tq(y)q(x')\0) . (11) 

Note that only the Fock terms appear since the Hartree (mean field) terms have already 
been sorted. Consequently the resulting equation of motion reads 

(if - mq - G((a) + Z7 5 7i<7r» + ^l°r 3 )G(x - x') 

2 (12) 

= 5\x - x') - E(x - y)G(y - x'), 

where S(x — y) stands for the non-local Fock selfenergy that depends on G(x — y), and an 
integration over y is understood. By a Fourier transformation and an isospin decomposition, 

A ik = A 5 ik + A 3 T ik + A - T ik + A+T ik^ 

t± = -j=(n±iT 2 ), (13) 



we obtain a Gor'kov [37| type equation, 

7°(w-/ i ± AiI /2) + S°±S 3 -G(7t) ?7 5 + v/2S^ W G° ± G 3 \ _ (l 
-G(tt)z7 5 + 7 (w-/iT/ii/2) + S TS 3 y y V2G ± ) ^0 

' (14) 

with h = a-k + (3(m q + G(a)) being the free single particle Hamiltonian with the constituent 
quark mass, M q = m q +G(cr). This form clearly indicates that the present subject is a pairing 
problem. The upper and lower double signs mean the u and d quark sector, respectively; 
both contain the same information. In the following we take the lower one. 

In order to solve Eq.( TH| and look into the spatial structure of the composite two body 



system, the pair wave function [38( given by the Bogoliubov amplitudes is necessary. The 
route is parallel to the non-relativistic case depicted in App. |A] This method was utilized 
for the nucleon pairing in Ref. [39]. In the present case, G° ± G 3 corresponds to the normal 
Green function and \/2G ± does to the anomalous one. First we express them in terms of 
the densities. The relativistic free quark field of 2-th flavor without pairing is expressed as 

= TTf EK^( ks ) e ~^' + KlV a (ks)e^} (15) 

v ks 



with k° = Ef. = ^k 2 + Mg. The number of single particle states must be doubled so as 
to have two energy states mixed by pairing interaction, as done by means of the Nambu 



representation [40| in field theoretical terms. The doubled states are diagonalized by means 
of the Bogoliubov transformation. Then the upper half states are regarded as unoccupied 
quasiparticle states while the lower half ones are occupied quasihole states. Therefore the 
particle states before transformation are regarded as superpositions of the quasiparticle with 
energy Sk and the quasihole with energy — Sk- Thus, in the present case, the quark field 
that defines G(x — x') is thought to be expanded in the same form as Eq.( fT5j) but with 

{+£k particle part with coefficient u % 
(16) 
—8k hole part with coefficient v\ 

with the Bogoliubov amplitudes specified below. Substituting it to Eq.(jS]) and Fourier 
transformation lead to 

G%{^ k) = £ U a (ks)U P M ( + { ff L) . ) 

' \cu — £k + ir] uj + S k — in/ 

(b% J>\ J , „ J> u s) 



+ J2v a (-k-s)Vp(-k-s)(- c . c . 

V uj — t k + ITj UJ + S k — If] 

+ > U a (\u)Vp{-k - S) [ . + 

\u — Sb + in 



(17) 



u) - Sk + in u + S k - if] 



+ ^7 a (-k-^( kS )( ^y + 

\u — Sb + in 



^ui — Sk + in u + Sk-ir]- 

Next, the densities such as (aa^) are expressed in terms of the Bogoliubov amplitudes by 
specifying the relevant transformation. In general, the exchange of quantum mesonic field 
produces non-local interactions of the type a l -a l , W-b 1 , tf-b 1 , a l -a\ b % -V , and a l -V (i ^ j). 
Therefore the quasiparticle takes the form of Eq. flB2l) . To be specific, however, here we 
consider a l —V that leads to the momentum dependent pionic gap function, and a l -a l and 
b l -b % that lead to the Fock mass, among them. Then the two types of quasiparticles specified 
in App. [B] decouple from each other; a 1 and b 1 ^ in q l become a constituent of different kind 
of quasiparticles. Then the normal and anomalous propagators are given as 

1 1 



s s 

V2G-p = ( E U a Vpu l v % * + E VjjpvW*) ( 



uj - Sk + in u + S k - in 
1 1 



uj — Sk + in ui + Sk — in 



(18) 



X 



0. 



Note that the expectation values arisen from the commutation relation are already sub- 
tracted in the backward terms. 

Substituting these expressions back to Eq. (IT4"l) and taking residues at uj = £ k , finally we 
obtain a 4 x 4 hermitian matrix equation at each k, 

( e-E k - fii/2 - m 2 -tt \ 

e + Ek — /ii/2 — m 2 —n 

-tt e + E k + fii/2 - mi 

-tt e - E k + m/2 -m x ) 

(A\ 

B 

C 

\ D J 

(19) 

Here the eigenenergy is denoted by e since both the quasiparticle and quasihole solutions 
are obtained from this, and use has been made of 

hU = E k U, hV = -E k V. (20) 

The real Bogoliubov amplitudes are defined as 

A = u 2 = (Ola^lO), B = v 2 = (0|&V|0>, 
C = -iv 1 = -i(Q\bl u rf\0), D = —iu 1 = —i(0\a u rf\0}, (21) 

and all quantities appearing in Eq. (jT9l) are real. Among them, 

7r(fc) = -iU(k)(-G(it)rf + y/2E + )V{k) 
fc(k) = -iV(k)(-G(it)rf + y/2E + )U(k) 

represent the momentum dependent pionic gap functions for the du and ud condensation, 
respectively, while 

m 2 (k) = -U(k)(E° -Z 3 )U(k), 
m 2 (k) = -V(k)CE° -Z 3 )V{k), 

(23) 

m x {k) = -V{k){Y? + Y?)V{k), 
mi (k) = -U{k)(Z° + Z 3 )U(k) 



(22) 
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do the Fock masses. The first term in each equation in Eq. (l22j) stems from the momentum 
independent pion condensation (it) of the meson system, which produces a strong momen- 
tum dependence, U^V = M q /E}-, and the second one from the non-local Fock selfenergy. 
This type of 4 x 4 matrix equation appears also in the cases of the relativistic 1 flavor pair- 



ing including the Dirac sea [39| and the non- relativistic 2 flavor pairing [41] . Since the Fock 
selfenergy at a momentum k is a function of A(k')-D(k'), the equations for all momenta are 
coupled. Actually, when evaluating each matrix element of S, a 4- momentum integration 
is necessary. For the energy integration among them, we make an instantaneous approxi- 



mation, that is, energy transfer— > as in previous works [19|, |20j, |39J. As for the remaining 
3-momentum integration, the BCS type calculation needs a cutoff in general. In the present 
case it is thought to be around the typical hadronic scale. Therefore we adopt that for the 
standard NJL model for simplicity. Solving the coupled equations selfconsistently deter- 
mines all the physical quantities: The Bogoliubov amplitudes, quasiparticle energies, and 
the mass and gap functions at each [x\. Then the pair wave functions and the coherence 
length are calculated from them. 

Now we proceed to numerical calculations. Parameters used are the current quark mass 
m q = 0.0055 GeV, the momentum cutoff A = 0.63 GeV, the pion decay constant f n = 0.093 
GeV, the pion mass = 0.138 GeV, the potential parameter in the linear sigma model A = 
4.5, and the quark-meson coupling G = 3.3. The momentum space < k < A is divided to 
100 equi-intervals for the coupled Newton method. Calculations are done for fi\ < where 
the du condensation dominates. The results depend on the parameters quantitatively but 
the qualitative behavior is robust; this will be confirmed later with respect to the behavior 
of the coherence length, which is of direct physical relevance. 

First, we check the meson masses under the present approximation in Fig. [TJ The cusp 
just after the transition = m n is brought about by the neglect of the Coriolis coupling 
term in Eq.(J7J). Definitely, the eigenvalues of the 2x2 diagonalization after that in the polar 
coordinate representation are 

M 2 = ^(2A 2 «a> 2 + (tt) 2 ) + £ ± y/[2X^-(^) + ^ + 16\*(a)^), (24) 

while two non-zero eigenvalues of the 3x3 diagonalization in the Cartesian coordinate 
representation 26j are 



M 2 = \{2\\(a) 2 + (vr) 2 ) + 5/x 2 ± y/[2\*((a)* - - 3/x?]' + 16\*(a)*(n)*). (25) 



The present result given by Eq. (l24p . the lower one of which tends to when (71") approaches 



0, is not consistent with the one obtained in the frame of the chiral perturbation [42|, but 
this difference is a trade-off for obtaining the meson mixing. Practically, its influence is 
limited to just after the transition. 



> 
o 




0.4 0.6 

abs(Hj) (GeV) 

FIG. 1: (Color online) Meson masses given by the linear sigma model with the approximation 
described in the text. Note that it and 9 correspond to ix\ and iT2, respectively, at \fi\\ < m n . 

Figure [2] shows the results at \fii\ = 0.5 GeV ^> m^. Figure[2](a) is the quasiparticle energy 
diagram as a function of the relative momentum k (dispersion relation). Its unperturbed 
structure is quite simple: The positive and negative energy u (d) quark levels with i-E^ 
are shifted upward (downward) by |/ii|/2. Then, the negative energy u, that is the hole 
state of the u, and the positive energy d interact around the Fermi surface. This means the 
du pairing. Hereafter we name these quasiparticle (hole) levels the first, second, third and 
fourth, from the bottom. The third level, the lower quasiparticle, is the main interest in the 
following discussion. This lower quasiparticle consists only of A and C. In the usual pairing 
problem, for example in the case of Ref. 



this type of 2 x 2 equation can be cast into the 
form of the gap equation. In the present case, however, n(k) is represented as a function of 
A(k') and C(k') as 

ir(k) -- 



~ J2(v(k, k')2A(k')C(k') + v'(k, k')(A 2 (k') - C 2 (k'))), 
w 

n(k') 



2A{k')C(k') 
A 2 (k')-C 2 (k') 



e(fc') 



E w + £ + 



mi(fc / )+m 2 (fc / ) ' 
2 

mi (k')—m,2 (k') 
2 



mi(fc')+m 2 (fc') 



(26) 



e(k>) - 

Therefore the v' term due to the a-ir mixing prevents one from casting Eq. ffT9]) into the 
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form of the gap equation. Nevertheless, the notion of the pair wave function 38j] is useful 
for looking into the physical contents since A 2 — C 2 is small around the Fermi surface. 
Figure [2](b) shows the Bogoliubov amplitudes A and C. Aside from the bump around k = 
mentioned below, the hole character changes gradually to the particle character around the 
Fermi surface as the usual Cooper pairing. This leads to the peak in the pair wave function 
4>{k) = A(k)C(k) (see Eq. (l26|) ) shown in Fig. [2] (c). The bump around k = is a novel 
feature of the present case; this is brought about by the mesonic contribution (n) to the gap 
function ir(k) (see Eq. (1221) ) as shown in Fig. [2] (d). In this gap function, the mesonic and 
the Cooper pair components are comparable around the Fermi surface, whereas the former 
is dominant around k = because of the k dependence oc M q /E k . 



> 
a 
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FIG. 2: (Color online) Momentum dependence of various quantities at fii = —0.5 GeV: (a) the 
quasiparticle energies, (b) the Bogoliubov amplitudes, (c) the pair wave function, and (d) the gap 
function. Note that (b) - (d) are associated with the third (from the bottom) solution in (a). 



Figure |3] shows the /ii dependence of various quantities. Figure [3] (a) shows the pair 
wave functions at several /j,is as functions of the momentum. This shows that, leaving room 
for possible error related to the discussion about Fig. [fl at low the peak due to the 
Cooper pairing can not be seen. Actually, q and q are bound to each other for \fii\ < 2M q 



11 



/ <f)(k)jo(kr)k 2 dk. (27) 
Jo 



as shown in Fig. [3] (b). Thus, we can conclude that the pionic condensation has a mixed 
character: Purely bosonic just after the appearance of the condensation, then the Cooper 
pairing gradually grows as |/ii| increases with retaining significant bosonic component. To 
look into the spatial structure of Cooper pairs more closely, we Fourier transform <p{h) as 

2vr 2 " 

The results for several /iis are shown in Fig. |3] (c) as functions of the relative distance. 
Obviously those for higher |/ii| wave till longer distance. Figure US (d) graphs the coherence 
length, 

and [3] (e) the gap at the Fermi surface as functions of fi\. The obtained coherence length at 
low | is consistent with the value obtained by an analysis of the ti-tt scattering, (r 2 )^ = 
0.61 ± 0.04 fm 2 43J. In relation to heavy ion collisions, this value is very close to the 



typical inter-pion distance d at the freeze-out: An example of numbers, the charged particle 
multiplicity N c = 555 |44j and the source size V = (6.48fm) 3 30], and the fact that the 
pion is the most abundant, lead to d ^ (V/A^) 1 / 3 = 0.79 fm. The picture of a gas of bound 
mesons may apply to £ < d while that of a liquid (see also Ref. 45|) of Cooper pairs would 
be appropriate for £ > d although the latter realizes at rather high |/ii|. Figure [3] clearly 
indicates that the Cooper pairing becomes weakly coupled as |/ii| increases. Comparing these 
figures with corresponding ones in Ref. 19| , one can see that the Cooper pairing part of the 
present case is more weakly coupled than the case of color superconductivity, as represented 
by the narrower peak in <p(k) and longer spatial extent. Figure |3] (d) also shows the cutoff 
dependence; the dependence is weak. 

At higher |/ij|, in the present calculation |/ii| > 0.8 GeV, a gapless pairing (e < 0) takes 
place. The gapless dispersion is known to occur in the case of pairing between particles with 
different masses J^. In the present case, the Fock term produces the difference in the mass 
(see the denominator in Eq. (I26j) ) . 

Finally we look into the character of the fourth level, the higher quasiparticle, that 
corresponds to the Dirac sea pairing in Ref. 39j|. This level is of almost pure u quark 
particle character (D(k) ~ 1) for k ^ 0.1 GeV; but the d component strongly mixes around 
k = because of two reasons: 1) (tt) equally contributes to 7c(k) and ir(k) (but with the 
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FIG. 3: (Color online) Isospin chemical potential dependence of various quantities: (a) the k 
space pair wave function, (b) the twice of the constituent quark mass, (c) the r space pair wave 
function, (d) the coherence length, and (e) the gap at the Fermi surface, (d) also contains the 
cutoff dependence. 

opposite sign), and 2) the unperturbed energy difference between u and d is the same as 
that between u and d at k = 0. 

To summarize, we have studied the momentum dependence of the pionic gap function 
ir(k) that determines the spatial structure of the condensation by adopting the linear sigma 
model as an inter-quark interaction at finite isospin chemical potential as a first step towards 
the study of the asymmetric matter in the real world. Although confinement is not taken 
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into account in the present study, the character of the condensation is bosonic at low \fii\, 
then the Cooper pairing gradually grows as \fi\\ increases. This q-q pairing is weaker than 
the q-q pairing of the case of color superconductivity. The spatial structure (wave function) 
of the composite pionic system is expected to be measured in lattice QCD simulations as well 
as the Hi dependence of the magnitude of the condensation as signatures of the BEC-BCS 
crossover. The spatial structure may affect the description of pions created in heavy ion 
collisions. 



Appendix A: The Gor'kov formalism 



Gor'kov [37] first proposed a field theoretical method to describe the pairing problem. In 
addition to the normal Green function G(x — x'), the anomalous Green function F'(x — x') 
of (Tfy'ip ')) type is introduced there. The equation of motion of their Fourier transforms 
is given by 

iA U + \ft(w,k) J \0J 
where and A are the single particle energy measured from the Fermi surface and the 
momentum independent pairing gap, respectively. Its solution is 



u-£ k ui + £ k 

t -iA/2£ k iA/2£ k 

Fl (u, k ) = 7T- + —7. , 

to - £ k u + £ k 

«K=(^) 2 , 4 = J% + A 2 . (A2) 



Substituting them back to Eq. flAll) gives 

[(u-£ k )u k -Av k ]u k _|_ [{uj-£ k )v k +Au k ]v k 
w-£ k u+£ k 

j [A%-K^K]ut 1 j [Av k +(uj+£ k )u k ]v k 
w-£ k tu+£k 

The residues at u — £ k (quasiparticle) lead to 



)-l 









(A3) 



&A h= e »ri- <A4) 

a J \v k j yv k 

Those at u = —£ k (quasihole) lead to the same equation. Therefore the equation for the 
Green functions and that for the Bogoliubov amplitudes are equivalent. 
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Appendix B: The Bogoliubov transformation 

Replacing the spin a =t / 4 and a_k in the non-relativistic pairing problem by the isospin 
u/d and fc_ k , respectively, we obtain two Bogoliubov transformations relevant to the present 
case, 




at each momentum and spin. 

Since there is i^ 5 !^ between two flavors, here we take u 1 and v 1 are imaginary, u 2 and v 2 
are real. Then the two types of quasiparticle, rjl and ry d , can be represented collectively as 

2 

i=i 

where i — 1/2 correspond to w/gL 
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